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Abstract 

Context. The synchrotron-self Compton (SSC) radiation process is widely held to provide a close representation of the double peaked 
spectral energy distributions from BL Lac Objects (BL Lacs), which are marked by non-thermal beamed radiations, highly variable 
on timescales of days or less. Their outbursts in the y rays relative to the optical/X rays might be surmised to be enhanced in BL Lacs 
as these photons are upscattered via the inverse Compton (IC) process. 

Aims. From the observed correlations among the spectral parameters during optical/X-ray variations we aim at predicting corre- 
sponding correlations in the y-ray band, and the actual relations between the y-ray and the X-ray variability consistent with the SSC 
emission process. 

Methods. Starting from the homogeneous single-zone SSC source model, with log-parabolic energies distributions of emitting elec- 
tron as required by the X-ray data of many sources, we find relations among spectral parameters of the IC radiation in both the 
Thomson (for Low energy BL Lacs) and the Klein-Nishina regimes (mainly for High energy BL Lacs); whence we compute how 
variability is driven by a smooth increase of key source parameters, primarily the root mean square electron energy. 
Results. In the Klein-Nishina regime the model predicts for HBLs lower inverse Compton fluxes relative to synchrotron, and milder 
y-ray relative to X-ray variations. Stronger y-ray flares observed in some HBLs like Mrk 501 are understood in terms of additional, 
smooth increases also of the emitting electron density. However, episodes of rapid flares as recently reported at TeV energies are 
beyond the reach of the single component SSC model with one dominant varying parameter. Furthermore, spectral correlations at 
variance with our predictions, as well as TeV emissions in LBL objects (like BL Lacertae itself) cannot be explained in terms of the 
simple HSZ SSC model, and in these cases the source may require additional electron populations in more elaborate structures like 
decelerated relativistic outflows or sub-jet scenarios. 

Conclusions. We provide a comprehensive benchmark to straightforwardly gauge the capabilities and effectiveness of the SSC radia- 
tion process. The single component SSC source model in the Thomson regime turns out to be adequate for many LBL sources. In the 
mild Klein-Nishina regime it covers HBL sources undergoing variations driven by smooth increase of a number of source parameters. 
However, the simple model meets its limits with the fast/strong flares recently reported for a few sources in the TeV range; these 
require sudden accelerations of emitting electrons in a second source component. 

Key words, galaxies: active - galaxies: jets - BL Lacertae objects: general - radiation mechanisms: non-thermal 



' 1. Introduction diation (synchrotron-self Compton, SSC; see, e.g., Marscher & 

GearQIMI Maraschi et al.[[9§2>. 

. . . Blazars are among the brightest Active Galactic Nuclei (AGNs) ^ nl , ,, , j • „ c 

. , . - . „ , . .. r The BL Lacs are usually classified in terms or the frequency 
observed, with inferred (isotropic) luminosities up to L - c ■ , . , ; n , . ■ 1 1 nnci c 
,_47 _i . , «^xr of their synchrotron peak (see Padovam & Giommi 1995): for 
10 erg s . Actually these sources are widely held to be AGNs , t DT , nm u , , ,. • • t , 
,. . & _ , , 3 . , , the low frequency BL Lacs (LBLs) the peak lies in the infrared- 
radiating from a relativistic let closely aligned to the observer .. , , j . „ . . , 

. . , ,. . J , , ,. • • , , optical bands, while the inverse Compton (IC) component peaks 

line of sight; this emits non-thermal radiations with observed f » . * , ■ ■ . , , c OT T /Tir> T •> f 

rr r , , at MeV energies; instead, high frequency BL Lacs (HBLs) fea- 

nuxes enhanced by relativistic effects so as to often overwhelm . c . , ■ v if j j »u j ., mr , 

3 ... ture a first peak in the X-ray band and the second one at 100 GeV 

other, thermal or reprocessed emissions. . , , . . . c nT T , TT)T s 

r energies or beyond; intermediate frequency BL Lacs (IBLs) 

BL Lacs in particular are Blazars showing no or just weak stan d in the middle, with the first peak at optical frequencies 

emission lines. Their spectra may be represented as a contin- anc j the second peak at ~ 1 GeV. 

uous spectral energy distribution (SED) S v = vF v marked by „ T T „, . , , . , .... , . 

, , DJ . ... . . 3 BL Lac Obiects also show rapid variability on timescales of 

two peaks: a lower frequency one, widely interpreted as syn- , , . , , . , iL j a 

, f . . . , 3 , , ■ , days or shorter, during which the sources undergo strong flux 

chrotron emission by highly relativistic electrons; and a higher . .. f ua „ ° ° 

r \ I i- j j. i „ ° ^ variations often named flares . 
frequency counterpart, believed to be inverse Compton upscat- 

tering by the same electrons of seed photons either emitted by The homogeneous single zone (HSZ) SSC model (based, 

an ex ternal source (external Compton, EC, see e.g. Sikora et al. in particular, on a single electron population) constitutes an at- 

[T994j l, or provided by the photons of the very synchrotron ra- tractively simple source structure worth to be extensively tested 

and pushed to its very limits, as we discuss below; following 
Katarzyhski et al. (2005), Tramacere et al. (2007), we focus 
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IC fluxes of BL Lacs to be expected from an updated version 
of SSC. We use a realistic log-parabolic energy distribution of 
the electrons radiating the synchrotron photons and upscattering 
them by IC in either the Thomson or the Klein Nishina (KN) 
regime. 

2. Background relations 

2. 1 . Log-parabolic electron distribution and spectra 

The SEDs of BL Lacs exhibit two distinct peaks or humps. To 
describe these in detail a number of spectral models have been 
used over the years: power-laws with an exponential cutoff, bro- 
ken power-laws, or spectra curved in the shape of a log-parabola. 
The latter are being increasingly used, as they generally pro- 
vide more accurate fits than broken power-laws in the optical 
X-ray bands, with residuals uniformly low throughout a wide 
energy range (see discussion by Landau et al. I1986I Tanihata et 
al. I20D41 M assaro et al. I2004al Tram acere et al. l2007l Nieppola 
et al. I20061 Donnarumma et al. EOOQlFl 

So we adopt log-parabolae in view of their optimal perfor- 
mance as a precise fitting tool. Moreover, they are also effective 
in yielding simple, analytical expressions (given in their exact 
form in AppendixlAl following Eqs. IA.6IIA.9l and rA.15l ) for peak 
frequencies and fluxes in both in the inverse Compton regimes, 
Thomson and Klein-Nishina. Finally, the physically interesting 
feature of log-parabolic spectra is their direct link to the accel- 
eration processes of the emitting particles; in fact, as shown in 
Appendix|B] log-parabolic emitted spectra simply relate to log- 
parabolic particle energy distributions, that in turn naturally arise 
from systematic and stochastic acceleration processes as dealt 
with by the Fokker-Planck equation. 

Thus for the synchrotron SED we write 



\ vo 



(a. t -l)-Mog(^) 



(1) 



here a s gives the constant component of the spectral index for 
the energy flux and b s is the spectral curvature. In the following 
we shall denote with £ the frequency where the SED peaks and 
by S the related flux. Spectral curvature values are observed to 
be around 0.2 (Massaro et al. |2004bl l and decreasing with time 
and/or frequency as given by Eqs. IB. 1 fl and IB. 13l of Appendix 
IbI since the decrease is mild during flares (which are essen- 
tially dominated by systematic acceleration processes), we can 
consider a nearly constant spectral curvature to a first approx- 
imation. Then one finds the scaling relations £ oc By 2 6 and 
S oc B 2 y 2 nR 3 5 4 , where n is the particle density, R is the emit- 
ting region radius, B is the magnetic field, 6 is the beaming fac- 
tor, and y p m c 2 is the square root of the mean quadratic energy 
(thereafter r.m.s. energy) of the particle distribution. 

An analogous expression holds for the inverse Compton ra- 
diation 



■{fic-Vj-b, 



(2) 



which in the Thomson regime the IC has less curved spectra 
with b c w b s /2. At nearly constant curvature the parameters to 



1 The log-parabolic particle energy distribution could be made asym- 
metric by adding a powerlaw tail at low energies (see Massaro et al. 



2006) to represent also low frequency data, e. g., radio in addition to 
optical/X-ray data; however, here we will focus on the SEDs around 
their peaks and so we do not use such additions. 




Figure 1. Two examples of exact, steady SEDs provided by nu- 
merical simulations (code by Massaro 2007) after the HSZ SSC 
model from relativistic electrons with a log -parabolic energy dis- 
tribution N(y) oc y- s - rX °&y (see AppendixlBt. The two character- 
istic humps correspond to synchrotron emission (left) and to in- 
verse Compton scattering (right), either in the Thomson regime 
(blue line) or in the KN regime (red line). Standard sources that 
are so fitted include the low state of 0716+714 (an IBL, blue 
crosses) observations of November 1996 by BeppoSAX, AIT and 
EGRET (Tagliaferri et al. 12003]) ; and 1ES 1553+113 (a HBL, 
red crosses) observations of April 2005 by Swift and April-May 
2005 by MAGIC (Massaro |20U71 ). 



focus on are e oc B y p 5 and C oc B <y p n 



A p 5 and C oc B 2 y 4 p n 2 R 4 6 4 for the SED peak 
frequency and flux, respectively. In the extreme Klein-Nishina 
regime we obtain instead more curved spectra with b c m 5 b s , 
and e cc y p S and C oc Bn 2 R 4 5 4 hold (see AppendixlAlfor more 
details, and Massaro et al. l2008al[200gbb . 

Such spectra are radiated by electron populations with a log- 
parabolic energy distribution, namely 



My) 



N U 



— s— J" log I 



(3) 



The result is easily inferred and straightforward to compute on 
using the delta-function approximations for the single particle 
emissions (see AppendixFAl; it is found to hold also for the exact 
shape of the latter with the numerical simulations outlined in Fig. 
[1] and reported in full detail in Figs. IB.ll and fB.2l 

2.2. Source parameters 

Using the previous scalings, with the relation R = c At 5/(1 + z) 
between the observed variation time Af and the size R of the 
emitting region (assuming spherical symmetry in the rest frame), 
we obtain relations of source parameters with the observables. In 
the Thomson regime we find 



7 P K 



R oc £ -- S-- Afl f -1 C - * 



(4) 
(5) 
(6) 
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6 oc S s^ 1 C~* Ar J (7) 

nccfCie-iS^-Ari. (8) 
For the extreme KN regime we find instead: 

y p oc^i ei C5 Afis - * (9) 

Z? cc e l S§ Af^-I c~5 (10) 

fioc£§S§e f C~5 AT* (11) 

tfoce'S'f^C"? Af' (12) 

11 7 4 9 1 

n oc £7 c$ e"s s -5 Ars. (13) 



Note that a delay of the order of Af may elapse between syn- 
chrotron and IC spectral variations. 

Numerical values for the prefactors are given in Appendix 
ICl where normalizations are chosen referring to a typical 
LBL source with S ~ 1CT 11 erg ctrT 2 s" 1 , £ ~ 10 l4 Hz, 
C ~ 10" 11 erg cm" 2 s" 1 , e ~ 10 22 Hz; and to a typical HBL 
source with S ~ 1(T 10 erg cm" 2 s _1 , £ ~ 10 18 Hz, C ~ 
10~" erg cm" 2 s _1 , e ~ 10 26 Hz. Values of physical parame- 
ters for some sources are given in Tab. Q] Note that the simple 
Eqs. © - ([8j provide parameter values in close agreement with 
the ones needed to fit the observations of LBL or IBL sources 
in Thomson regime; on the other hand the extreme KN limit 
underlying Eqs. <j9j - ( [T3l would yield disagreement for HBLs, 
indicating that these sources just border the KN regime. 

We evaluate the position f-j of the synchrotron peak when 
the transition between the two IC regimes occurs, to find 

^^• 1510l 1bTGF(^) (1+zrlHz ' (14) 

or equivalently 

- 1.96 10" (^j" 1 10K>-i)(A) ( i +z) -i Hz. (15) 

So we conclude that in LBLs, which feature the synchrotron 
peak at optical/IR frequencies with £ w 10 14 Hz, the majority of 
photons of the IC peak are upscattered in the Thomson regime; 
whereas in extreme HBLs, which show the synchrotron peak in 
the X-ray band with £ as 10 18 Hz, the majority of IC peak pho- 
tons are upscattered in the Klein-Nishina regime. 

2.3. Spectral correlations during X-ray spectral variations 

Correlations in the synchrotron emission between S and b s with 
£ can be used to pinpoint the main driver of the spectral changes 
during X-ray variations (see also Tramacere et al. 2007). In fact, 
the synchrotron SED peak scales (for nearly constant spectral 
curvature as anticipated in 92.lt as S ocn R 3 y 2 B 2 6 4 at the peak 
frequency itself scaling as £ oc y 2 B 6. 

Thus on writing the dependence of S on f in the form of a 
powerlaw S oc we expect ar = 1 to apply when the spectral 
changes are driven mainly by variations of the electrons r.m.s. 
energy; a — 2 for dominant changes of the magnetic field; a = 4 
if changes in the beaming factor dominate; a = oo formally ap- 
plies for changes only in the number density of the emitting par- 
ticles. Results by the latter authors focus on y p and B as domi- 
nant drivers. Starting from such correlations we aim at predicting 
the expected correlations and flux variations in the y-ray band. 
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Particle r.m.s. energy variations are of particular interest be- 
cause they naturally arise as a consequence of systematic plus 
stochastic acceleration of the electron population (details in 
Appendix|B]). 

3. y-ray spectra 

3.1. Spectral correlations 

The above results lead us to expect specific correlations between 
SED peak frequencies and fluxes as follows. In the case of dom- 
inant r.m.s. particle energy variations we expect for the syn- 
chrotron emission 

Sccf, (16) 
while for the IC we have 

Coce (17) 
for LBL objects, and 

C ~ const (18) 

for extreme HBLs. 

In the case of magnetic field variations we have for the syn- 
chrotron emission 

Soc£ 2 , (19) 
while for the IC we have 

Coc £ 2 (20) 
for LBL objects, and 

C oc £ T (21) 

formally with r — oo for extreme HBLs. 

To complete the above picture, we consider how saturation 
affects correlations. Denoting the total number of particles with 
Nt, in the jet frame we have for the synchrotron luminosity 

L s oc Nt y 2 ,, while for IC luminosity we have Lq oc (Nt y 2 ,] and 

Lq oc Nj in Thomson and KN regime, respectively. Consider 
the case where the total power available to the jet is limited, e. 
g., in the BZ mechanism (Blandfor & Znajek [T977t of power 
extraction from a maximally rotating black hole with an ex- 
treme, radiation pressure dominated disk that holds the pressure 
of the magnetic field facing the black hole horizon (see discus- 
sion by Cavaliere & D , Elia |2002| and references therein). Then 
the power has an upper bound Lbz ~ 2 10 45 M9 erg s _1 ; as a con- 
sequence, we expect an effect of saturation on the SED peaks 
as the emitted power, which represents a substantial fraction of 
the total jet power (see Celotti & Ghisellini 2008), approaches 
Lbz- For the total power, including radiative, kinetic and mag- 
netic components, L m ~ Lbz we envisage Nj ~ Jj 2 , that is, the 
particles effectively accelerated to the increasing y p decrease in 
number on approaching the BZ limit. 

The data analyses by Tramacere et al. (2007 2009) for the 
variations of Mrk 421 are consistent with the beginning of such 
a saturation effect, as they show a declining value of a when 
the peak moves toward the highest energies and approaches 
L a 10 46 erg s _1 . Partial saturation effects can also be explained 
in terms of systematic and stochastic acceleration processes; in 
particular, when the stochastic prevails over the systematic ac- 
celeration (as is bound to occur at high energies, see Appendix 
lAt we expect a to decrease toward 0.6, i.e. S oc £°- 6 . 
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Source name 


z 


7p 


/J (cm) 


B(G) 


S 


n(cm i ) 


r 


0716+714 (low state) 


0.31 


3.36 10 J 


2.61 10" 


1.51 10" 1 


10 


0.41 


1.18 


0716+714 (high state) 




9.92 10 3 


7.79 10 17 


7.64 10~ 2 


10 


0.03 


1.18 


1ES 1553+113 


0.25 


3.39 10 4 


0.64 10 1S 


9.80 10-' 


10 


2.04 


1.49 


Mrk501(low state) 


0.034 


1.88 10 5 


0.38 10 16 


9.29 10~ 2 


15 


0.62 


0.81 


Mrk 501(high state) 




2.58 10 5 


0.59 10 15 


3.85 10-' 


15 


23.49 


0.74 



Table 1. Source parameters for some BL Lac sources. 



3.2. Related variabilities 

We now focus on related variabilities between IC and syn- 
chrotron fluxes; here we are mainly interested on increases of 
peak frequencies and fluxes, as driven by acceleration processes 
of the emitting particles or by a growing magnetic field (see also 
Katarzynski et al. 2005); in Appendix |Pl we briefly discuss the 
problems posed by the decay phase as stressed by Katarzynski 
et al. 

For LBLs, that radiate mostly in the Thomson regime, we 
predict from Eqs. ( fToT l - (f2Tb a quadratic or linear variation of C 
with respect to S , depending on the parameter that mainly drives 
source variations. That is to say, we expect 

AlnC~2AlnS (22) 

in the case of dominant particle r.m.s. energy variations, or 

AlnC-AlnS (23) 

for dominant magnetic field variations (see also Figs.|2]and|4|. 

Instead, for HBLs that radiate closer to the Klein-Nishina 
regime, we obtain a weaker y-ray variability due to the decreas- 
ing KN cross section, leading in the extreme to 

AlnC~0 (24) 

in case of dominant particle r.m.s. energy variations, or 

AlnC~-AlnS (25) 
2 

for dominant magnetic field variations. Expected correlations are 
summarized in Table [2] 

From Figs. [2] and [4] it is seen that the transition to the KN 
regime has three effects on the IC spectrum: first, it brakes the 
peak frequency increase, because the energy the photon can gain 
is limited to the total electron energy ymc 2 ; second, it reduces 
the flux increase at the SED peak, as a consequence of reduced 
cross section; third, it increases the spectral curvature near the 
peak as a consequence of frequency compression. 

Viceversa, variations observed in IC section of the spectrum 
are related via Eqs. d22b - d25l ) to variations expected in the syn- 
chrotron emission; in particular, in HBLs y-ray variations ought 
to have enhanced counterparts in X rays. 

Up to now Figs.|2]and|4]may be interpreted as a collection of 
different sources or of unrelated states of a single source; hence- 
forth we will focus on the interpretation in terms of evolution in 
a flaring source where the parameters vary in a continuous fash- 
ion; in particular, y p increases under the drive of systematic and 
stochastic electron accelerations, in keeping with the descrip- 
tion in terms of a continuity or kinetic equation of the Fokker- 
Planck type as given in detail in AppendixlBl Signatures of such 
an evolution are provided not only by continuous growth of the 
peak frequencies, but even more definitely by the irreversible de- 
crease of the spectral curvature (in terms of time or frequency) 
under the drive of the stochastic acceleration, see Eq. ( IB- 1 3D and 




log v (Hz) 



Figure 2. To illustrate the differences of the IC SEDs in the 
Thomson and in the extreme KN regime, we show the results 
of numerical simulations (code by Massaro 2007) after the HSZ 
SSC model, for different sources with larger and larger r.m.s. 
energy of the radiating particles and other parameters kept con- 
stant (including the curvature of the particle energy distribution). 
The vertical line indicates the frequency of the synchrotron peak 
where the IC scattering changes over the Thomson to the KN 
regime (see Eq. ( fT4b ). Actual sources in their evolution (repre- 
sented in Fig. IB. II and IB. 21 ) vary their spectral curvature some- 
what and only span a limited range in frequency and flux to the 
left (LBLs) or to the right (HBLs) of the line, with only the IBLs 
likely to cross it. 



the observations by Tramacere et al. ( 120071 [2009b . In fact, for a 
single electron population the curvature is to slowly decrease or 
stay nearly constant on the timescale of the systematic accelera- 
tion. So a sudden increase of curvature signals the injection of a 
new population/component. 

In the following we will neglect radiative cooling, which 
does not strongly affect the spectral shape around the peaks; such 
a process will be particularly relevant when the source, after flar- 
ing up, relaxes back to a lower state upon radiating away its ex- 
cess energies from high frequencies downwards (see Appendix 

m. 

3.3. Specific sources 

We have applied the above simple expectations to two sources, 
0716+714 (IBL) and Mrk 501 (HBL), among the few to date to 
provide extended simultaneous coverage of two different source 
states in X rays and y rays. 



A. Paggi et al.: SSC radiation in BL Lac sources, the end of the tether 



5 



process 


peak flux and frequency 


flux-frequency correlation 


5 - C correlation 


synchrotron 


S ocR 3 B 2 y 2 ,n6 4 
? x By 2 p S 


S °z£"< 


a= 1 \y p ) 
, a = 2 (B) 




IC Thomson 


C oc R 4 B 2 y 4 n 2 S 4 
e°cBy 4 p S 


Coc e». 


r « = 1 (r P ) 
[ a = 2 (B) 


Coc5 2 
CocS 


ICKN 


CocR 4 Bn 2 S 4 
(Eocy p S 




r a = (y p ) 
a = oo (B) 


C ~ const 
CocS 112 



Table 2. Spectral correlations (for r s= const). We have denoted with (y p ) or (B) the variations driven by increases of rms electron 
energy or magnetic field, respectively. 



Alog(e)/Alog© 
Alog(C)/Alog(S) 
b/b. 

Mrk501 4/1997 
0716+714 11/1996-9/2007 



16 17 

log \ (Hz) 




20 22 
log(v) (Hz) 



Figure 3. Computed changes of spectral parameters correspond- 
ing to extended increase of r.m.s. particle energy with the 
other source parameters kept constant. As functions of £ — 
£(5/0.1 G) -1 ^(<y/10) -1 (l + z) (where £ is the synchrotron SED 
peak frequency, see Eq. (fT4li). the full line represents the ratio 
Alne/Aln£, the dashed line the ratio A In C/ A In S and the dot- 
dashed line the curvature ratio bjb c at the peaks. Circles rep- 
resent the parameters of 0716+714 as observed on November 

1996 and September 2007 (Tagliaferri et al. 120031 Chen et al. 
2008 ), and diamonds represent Mrk 501 observations of April 

1997 (Djannati-Atai et al. 1 19991 Massaro et al. |2006l l, each with 
the same colours as the lines. Changes of n only would yield 
A In C = 2 A In S at all frequencies, represented by the horizon- 
tal dotted line; added to a pure r.m.s. energy increase, a moderate 
increase of n can account for the upward deviations of observed 
values. 



The first is widely considered to be an IBL. In the low state 
observed by Automatic Imaging Telescope (AIT) in November 
14 1996, the synchrotron peak at about at about 10 14 5 Hz falls 
quite below the threshold frequency expressed by Eqs. ( TBI ) or 
(JT3J, so we expect for this source the IC scattering to occur 
mostly in the Thomson regime; simultaneous y-ray observations 
were carried out by EGRET. The high activity state of the source 
is described by data relative to GASP project of the WEBT and 
AGILE-GRID observations of 2007 September 7-12 (see Fig.|6>. 
The source has shown similarly low flux levels during EGRET 
observations 2003 ; so for flares with short duty cycle, it is not un- 
reasonable to evaluate variations between the two distant states 
in the absence of closer observations; curvature variation is not 
required to fit the data, as we expect at lower energies where sys- 



Figure4. As in Fig. [2] we show SEDs obtained from numer- 
ical simulations (code by Massaro 120071 ) of HSZ SSC radia- 
tions, for different sources with larger and larger magnetic field 
B and other parameters kept constant. The vertical line again in- 
dicates the transition frequency from the Thomson to the KN 
regime (see Eq. (fTBIl); as B increases, the IC scattering changes 
its regime, and reduced cross section yields similar compression 
effects as in Fig. [2] 

tematic dominates the stochastic acceleration (see AppendixlAl. 
As shown in Figs. [3] and [5] peak flux and frequency variations 
of 0716+714 appear to be in between the lines representing in- 
creasing r.m.s. particle energy and increasing magnetic field, and 
so they may be described in terms of HSZ SSC with simultane- 
ous variations of these two parameters. 

For the HBL source Mrk 501 we consider the two states de- 
scribed by Massaro et al. ( 120061 ), with simultaneous BeppoSAX 
and CAT observations of April 7 and 16 1997. To understand the 
behavior of this source, from Figs. [3] and [5] it is seen that vari- 
ations of B alone are ineffective and require complementary n 
increase by a factor of about 30; variations of y p are adequate 
in the KN regime, and only require a complementary doubling 
of n (which by itself would yield A In C w 2 A In S in both scat- 
tering regimes) to yield higher y-ray flux increase. Note that in 
going to higher energies the spectral curvature decreases some- 
what (from b = 0.161 ± 0.007 to b = 0.148 ± 0.005 for the low 
and high state, respectively) in keeping with model predictions; 
such a behavior is interesting as it marks a smooth growth, if any- 
thing, in the number of emitting particles as contrasted with sud- 
den, substantial re-injection of nearly monoenergetic electrons 
that would suddenly reverse the otherwise irreversible decrease 



6 



A. Paggi et al.: SSC radiation in BL Lac sources, the end of the tether 





Alog(e)/Alog© 




Alog(C)/Alog(S) 




b A 


♦ 


Mrk501 4/1997 


• 


0716+714 11/1996-9/2007 




log | (Hz) 




log v (Hz) 



Figure 5. Computed changes of spectral parameters as driven 
by an extended increase of the magnetic field with the other 
source parameters kept constant. The black full line represents 
the value of the ratio A In e/ A In the red dashed line rep- 
resents the value of the ratio A In C/ A In S and the blue dot- 
dashed line represents the ratio b s /b c , for different values of 
| = ^( 7/ ,/10 4 )10 (1 " 1/56)/2 (5/10r 1 (l + z) (where £ is the syn- 
chrotron SED peak frequency, see Eq. (fT5]l), data points are the 
same as in Fig. [3j compared with the latter considerably larger 
changes of n are needed to recover the observed values. 



of r and b. Such a smooth growth may obtains in a scenario of 
an expanding blast-wave that progressively involves more elec- 
trons (see, e. g., Ostriker & McKee [1988l Lapi et al. 120051 Vietri 
120061 . 

We stress that in HBL sources even small flux variations in y 
rays are expected to have enhanced and observable counterparts 
in X rays. We suggest that such variations should be checked 
upon y-ray "alarms", the inverse triggering relative to usual. The 
absence of such lower energies counterparts will indicate, for 
example, a flare driven dominantly by a particle number density 
increase associated with a magnetic field decrease n oc B 2 , caus- 
ing S ~ const but a decreasing synchrotron peak frequency; such 
a kind of flare may easily drown into the primary synchrotron 
component. 



3.4. Limiting timescales 

As shown in Appendix|B] systematic and stochastic acceleration 
processes occur in the sources on timescales t\ = and ti = 
1/A.2, respectively. Limits to the variability timescales Af relate 
to the apparent size of the emission region by Af > R (1 + z)/c 6 
(see jj2.2| i; rapid variations require small emitting sources, but 
these could be optically thick to the pair production process. In 
fact, y-ray photons collide with less energetic ones to produce e ± 
pairs, and the cross section of this process tops at a value around 
<Tr/5 (where crj is the Thomson cross section) when the y-ray 
photon frequency v and the target photons frequency v, satisfy 
the relation v, * m 2 c 4 6 2 /[h 2 (1 + zf v]. 

For the y rays to escape from the emitting region the source 
"compactness" intervenes and can set lower bounds on the 
beaming factor (Cavaliere & Morrison |1980t Massaro 120071 



Figure 6. Spectral fit of two different states of the IBL source 
0716+714. The low state (blue line) data are relative to 1996 
November 14 simultaneous observation by AIT, BeppoSAX and 
EGRET; the high state (red line) data are relative to 2007 
September 7-12 simultaneous observations by GASP-WEBT 
and AGILE-GRID (Tagliaferri et al. [20031 Chen et al. |2008l . 



Begelman, Fabian & Rees 2008). In fact the optical depth for this 
process is expressed in terms of observational quantities (primed 
quantities refer to the jet rest frame) as 



/ , \ cr T , , 9 err h D 2 a S v I v,\ 

rrr M = * f n ph W ) „ --^ ^ (1 + z) 4 - [ j) , (26) 

where D is the source distance expressed in Gpc, with the short- 
hand B = b log ( j J. Escaping y-ray radiation requires the source 

to be optically thin, that is, r yr \Y y J < 1, resulting in 



5> 



9 <r T h l+2 P 

20 m 2+2/3 c 6+4/? 



D 2 (1+z) 4 +2/J 



Ar 



(27) 



Thus rapid flux variations set a lower bound on the beaming 
factor; the highest bound between the previous relation and 
6 > R (1 + z)/c At is to be relevant. 

The exponent B contains a weak logarithmic dependence 
on 6 that may be neglected to a first approximation, and has 
a magnitude which depends on the source spectral properties; 
for definiteness, we use for b the typical value 0.2, and for 
v * 10 27 Hz. For a LBL with £ ~ 10 14 Hz we obtain B * 0.4, 
while for an extreme HBL with £ ~ 10 18 Hz, we obtain B * 
-0.5. A variation timescale Af ~ 5 minutes for a LBL with 
S ~ 10~ n erg crrT 2 s~' implies S > 30, while in an extreme 



HBL with S ~ 10~ 1U erg cm' 2 s" 1 we find 5 > 15. 

In the specific case of PKS 2 155-304 discussed by Begelman 
et al. (120081 1, the emitted power L =s 10 46 erg s~' corresponds to 
S ~ 4 10~ n erg cm" 2 s"\ with £ ~ 10 16 Hz, to give B ~ 0.1 
and 5 > 50 for Af ~ 5 minutes. We note that the source, even if 
widely considered a HBL, shows peculiar behaviors with respect 
to other TeV HBLs (see Massaro et al. 2008b); moreover, it does 
not satisfy the condition £ » £7- and so it cannot be considered 
an extreme HBL, as discussed below. 

Extreme values of the beaming factor 5 ~ 50 have been pro- 
posed to account for large peak separations (Konopelko et al. 
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2003 ) or to explain very rapid spectral variations (Begelman, 
Fabian & Rees 2008 ), formally with At ~ 5 minutes for the y-ray 
flares of PKS 2155-304 in 2006 July 28 (Aharonian et al. 12007) 
and Mrk 501 in June 30 and July 9 2007 (Albert et al. 12007a) . 
Note that with BL Lac spectra realistically curved, the observed 
flux variations may be enhanced due to a slope effect best legible 
on the differential flux F v \ that is, when fluxes are measured at 
frequencies where the spectral slope is steep (as may be the case 
for PKS 2155-304), a strong observed flux variation implies only 
a mild variation of the peak flux and requires smaller values of 

6. Otherwise, a peak variation on a scale At ~ 5 minutes would 

_ i 

require in the SSC model 6 ~ 100 for LBLs where 6 oc At '- (see 
Eq. (0), and 5 ~ 20 for an extreme HBL where 5 oc At~* (see 
Eq. ([13). 

As stated under 33.31 in HBL sources even small flux vari- 
ations in y rays should have enhanced counterparts in X rays. 
These may become hard to observe for example in a flare dom- 
inantly driven by particle number density increase associated 
with a magnetic field decrease (i.g., with n oc B 2 ), leading to 
a synchrotron emission easily drowned into other components. 
On the other hand, components with 6 > 20 appearing in an 
HBL spectrum go beyond the HSZ SSC model with variations 
of one dominant parameter, and therefore require a more elabo- 
rate source structure. 

In such cases the next natural scenarios are provided by 
decelerated relativistic outflows (Georganopoulos & Kazanas 
2003 ), or by nested spine-layer jets (Tavecchio & Ghisellini 
12008) and jets in a jet (Giannios et al. 12009) . 

4. Beyond the SSC 

In Table [2] we have provided a comprehensive benchmark to 
gauge the performance of the HSZ SSC model for flaring BL Lac 
objects. This allows us to easily recognize events that may be 
accounted for within the model with physical variations of key 
source parameters, or that instead require more complex source 
structures. 

Toward that purpose, we have used realistic log-parabolic 
spectral shapes produced by log-parabolic energy distributions 
of the emitting electrons, and studied IC radiation both in the 
Thomson and the Klein-Nishina regimes. 

In the model we expect S to dominate C fluxes; in fact, mov- 
ing to higher energies in a collection of different sources or in a 
prolonged evolution of a given source, we expect the emission 
to drift out of the pure Thomson and approach the KN regime, 
where the cross section decreases and limits the IC fluxes en- 
forcing C < S . Such a relation is found to hold in a number of 
sources, and in particular in HBLs (see for example Tagliaferri 
et al. l2U08l . 

From Table [2] we stress here the source spectral variations 
predicted in y rays. For LBLs, where the IC scattering mostly oc- 
curs in the Thomson regime, we recover the standard quadratic 
increase in y-ray fluxes with respect to IR-optical ones, ex- 
pressed by A In C ~ 2 A In S for dominant particle r.m.s. energy 
variations (see Eq.l22b. Instead for HBLs, in which the IC scat- 
tering approaches the Klein-Nishina regime, we expect smaller 
or even vanishing increases of the y-ray flux relative to X rays, 
that is, A In C - for dominant particle r.m.s. energy variations 
(seeEq.EU). 

Comparing with specific sources, we find the HSZ SSC 
model to be adequate for most LBLs at the present observa- 
tional levels; whilst for example the HBL source Mrk 501 in 
April 1997 showed a y-ray flux increase appreciably stronger 



than expected if it were dominated by just one driving parame- 
ter. We explain this behavior within the model in terms of an ad- 
ditional, smooth and moderate increase (by a factor around 2) in 
number density of the electrons responsible for the emission (see 
33.3) . Such conditions can still be provided by a single electron 
population; this should be marked by a continuously decreas- 
ing spectral curvature b, as indicated by current data (Massaro 
et al. 12006) . a feature providing a potentially powerful signa- 
ture to closely check on further data. Instead, for dominant mag- 
netic field variations we would expect A In C ~ A In S for LBL 
sources (see Eq. d23l . and A In C ~ j A In S for HBL sources (see 
Eq.l25); but in that case the required increase in particle number 
density should be much higher by a factor of about 30, hard to 
interpret in terms of a single electron population. 

On the other hand, on using the previous relations inversely, 
we see that even small y-ray variations in HBL sources ought to 
have enhanced counterparts in the X rays, unless flare activity 
were driven by an additional jet component, for example one 
with higher magnetic field but lower particle density; thus the 
corresponding emission is easily overwhelmed by, or drowned 
into the main synchrotron, but then the rest frame acceleration 
times must be short enough to involve a significant fraction of 
the emitting region. We suggest this is as a critical test for the 
simple model. 

Another limitation arises from rapid peak flux variations re- 
quiring large values of 6 as shown in 33.41 In the extreme, the 
few sources with particularly fast y-ray increases observed so 
far, like PKS 2155-304 in 2006 July 28 and Mrk 501 in June 30 
and July 9 2007 (see 33.4) , require additional components with 
very high beaming factors that go definitely beyond the simple 
SSC model. 

Finally, another problem for the model arises from LBL ob- 
jects showing substantial emissions in the GeV-Tev range that 
cannot be explained it terms of the simple HSZ SSC model 
(considering that second order IC scattering is negligible, see 
Massaro 2007) because the IC peak falls for these sources 
around 0.1 GeV; this may be the case for BL Lacertae itself 
(Albert et al. 12007b) and possibly similar sources like M87 
(Georganopoulos et al. 2005 ) and Cen A (Lenian et al. 2008; 
Aharonian et al. 12009) . 

In all these cases the source may require more elaborate 
structure, like decelerated relativistic outflows or sub-jet sce- 
narios (see Georganopoulos & Kazanas [2003; Tavecchio & 
Ghisellini 120081 Giannios et al. 2009). We stress that the in- 
jection of a second, nearly monoenergetic electron population 
is expected to be marked by a sudden increase of the spectral 
curvature. 



5. Discussion and conclusions 

We propose that the sources of BL Lac type may be conveniently 
ordered in a succession spanning from smooth variations of one 
or a few dominant SSC parameter, up to the appearance of truly 
different components that ultimately break through the limits of 
the simple HSZ SSC model. Even more so at increasing ener- 
gies and frequencies, that would imply weaker and weaker y-ray 
fluxes relative to X-ray, owing to KN cross section effects. This 
picture is supported by the following two lines of evidence. 

First, a recent statistical study of Third EGRET Catalogue 
data (Mukherjee et al. 120011 Casula 120081 Vagnetti, Casula & 
Paggi in prep.) shows for BL Lac objects observed at 30 MeV -f 
30 GeV a weak y-ray variability on average, compared to the 
FSRQ Blazars. Also the time-structure functions for these two 
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classes of objects indicate a similar trend, within the limitations 
of the sample. 

Second, the first AGILE-GRID Catalogue of high confidence 
gamma-ray sources (Pittori et al. 2009) apparently shows few 
if any new BL Lac sources other than those already detected 
by EGRET. This circumstance suggests the sources observed by 
AGILE to be more powerful than the average; it leads to expect 
for the majority of the sources either reduced flare activity in this 
band, or weak average fluxes. We relate these features to reduced 
Klein-Nishina cross section that tends to limit both the average 
fluxes and the flares. 

With Fermi, the first three months of observations (Abdo 
et al. 120091 ) have yielded a substantial number of new BL Lac 
sources, and in particular a larger fraction of HBLs close to the 
ratio known from radio surveys (Nieppola et al. 2006). This is 
consistent with our expectations, considering the better instru- 
mental sensitivity, especially at high energies generally favor- 
able to the HBLs; for these sources we suggest to test the model's 
outcomes: C lower than S , and generally small (though possibly 
fast) flux variations compared to the LBLs 

We plan to check our theoretical predictions against the pub- 
lic availability of the first year data from Fermi observations, 
and with simultaneous, multi-wavelength observations compar- 
ing the two basic subclasses of LBL and HBL of the BL Lac 
sources. 
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Appendix A: Log-parabolic spectra 

In this Appendix we show how electron populations with a log- 
parabolic energy distribution of the form expressed by Eq. (0, 
that is, 



N(y) 



-s-rlos — 
B \ro/ 



(A.l) 



emit log-parabolic spectra via the SSC process. The related par- 
ticle synchrotron emissivity 



/HP 
dyN(y) 



(A.2) 



is easily computed on using the close approximation to the 
single particle emission in the shape of a delta-function (see. 
Rybicki & Lightmann ll979l) . that is, ~ P s 5 (y - y 2 v c ) with 

P s = 1 /6tt ct t y 2 c B 2 and v e * 1.22 10 6 B Hz (with B measured 
in Gauss) is the synchrotron critical frequency. This leads to 
(Massaro et al. 2004a ) a log-parabolic differential flux 



Fi*F \-\ 
and to a SED again of log-parabolic shape 



v_<ft-i)-*,tag(£) 



(A.3) 



(A.4) 



its slope at the synchrotron reference frequency vo is given in 
terms of s, by 



s- 1 



the spectral curvature by 



(A.5) 



(A.6) 



and the peak value 5" oc p} B 2 ny 2 p y[r occurs at a frequency 

oc By 2 10" (y p , n, B and R are defined in ^2. II of the main 
text). 

For IC radiation in the Thomson regime we may write to 
a fair approximation ^ * p c S (y - | y 2 (see Rybicky & 
Lightmann ll9791 ) where P c = ^crry 2 c e v is the power radiated 
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by a single-particle IC scattering in the Thomson regime, having 
denoted with e v oc R B 2 n y 2 the synchrotron radiation densityQ . 
We obtain once again a log-parabolic SED 



Vvo 



-(a c -\)-b c log(v/vo) 



(A.7) 



where the slope at the IC reference frequency t>o is given by 
s-l 



and the spectral curvature reads 



r 



(A.8) 



(A.9) 



The peak value C oc R 4 B 2 n 2 y 4 yfr is attained at a frequency 
e' oc By 4 p . 

For the Klein-Nishina (KN) regime instead it necessary to 
consider the convolution 



f v =h j dyN(y) J dvvN,Kiv,r. y) 



(A. 10) 



where v and v are the electron frequencies before and after the 
scattering, respectively, N v is the number spectrum of seed pho- 
tons, and K(v, v,y) is the full Compton kernel (Jones 1968). Only 
in the extreme KN regime one may again approximate 



1 / hv , 

K{v,v,y)*—6\\ j ; 

y z \ ymc L ' 



(A.ll) 



on approximating N v with its mean value for a homogeneous, 
spherical optically thin source 



Ny * (Ny) 



3rj; 



4 c hv 

one obtains again a SED with a log-parabolic shape 

( \-(ae-l)-Mog(v/vo) 
ib) 

Now the slope at vo 

a c w s, 

is steeper, and the spectral curvature 

b c a r 



(A. 12) 



(A.13) 



(A. 14) 



(A.15) 



is larger than in the Thomson regime. The peak value C oc 
R 4 Bn 2 yfr occurs at a frequency e' oc y p \QT~r . 

The transition between the two regimes occurs when 
2 y ma x h ~ mc 2 , where y max — yjg' /v c , that is, when 

^^• 151 ° 15 (oTGH^) (1+zrlHz (A ' 16) 

holds, or equivalently 

£r * 1.9610 16 fef loiO-^^jd Hz. (A.17) 

2 Note that in Thomson regime no direct relation appears between 
electrons emitting peak synchrotron photons and those responsible of 
IC peak radiation (see Tavecchio et al. 1 1998b . 



We close these calculations with two remarks. First, the 
(primed) quantities used here refer to the rest frame of the emit- 
ting region, while the observed (unprimed) quantities must be 
multiplied by powers of the beaming factor 



1 



T(l -B cos6) 



(A.18) 



r being the bulk Lorentz factor of the relativistic electron flow 
in the emitting region, 9 the angle between its velocity v and the 
line of sight, and B — v/c; so we obtain S = S' 6 4 ,C — C 6 4 , and 
% = 5, € — e' S. Second, we note that the actual specific powers 
4^ radiated by a single particle differ from delta-function shape, 
slightly for the synchrotron emission and considerably for the 
IC radiation in the Thomson regime (see discussion in Rybicky 
& Lightmann 1979). However, the above spectral shapes still 
approximatively apply, as confirmed by numerical simulations 
(Massaro 2007); in detail, the convolution with a broader single 
particle power yields a less curved spectrum, so that equations 
IA.6I and [A~9l become b s m r/5 (Massaro et al. 2006) and b c =s 
r/10, respectively. 



Appendix B: Particle acceleration processes 

Here we derive a log-parabolic electron energy distribution 
Af (y, t) from a kinetic continuity equation of the Fokker-Planck 
type; following Kardashev ( 1962), in the jet rest frame this reads 



dN 
~dt 



d 

-Ai(t) — (yN)- 
ay 



dN 
dy 



(B.l) 



where ym c 2 is the particle energy, t denotes time, and Ai and A2 
describe systematic and stochastic acceleration rates, occurring 
on timescales t\ = l//li and ti = IM2, respectively. For exam- 
ple, in the picture of Fermi accelerations (e.g., Vietri 2006), the 
accelerations rates A\ and Ai can be expressed in terms of phys- 
ical quantities related to processes occurring in shocks; in this 
framework a plane shock front of thickness l s moves with speed 
V s and gas clouds of average size £ move downstream of the 
shock with speed V; it is found that A\ = V s /{ s and Aj - V 2 /2c{ 
hold (see Kaplan [T9561 Paggi l2007l l. A 2 <k Ai. On the other 
hand, numerical values of A\ and A2 are directly derived from 
the emitted spectrum as shown later. 

The Fokker-Planck equation ( IB.lt describes the evolution 
of the electron distribution function; with an initially mono- 
energetic distribution in the form of a delta-function N(y, 0) = 
n 6(y - yo) (n is the initial particle number density) the solution 
at subsequent times t takes the form of the log-parabolic energy 
distribution assumed in Eq. © of the main text (see also IA.U . 
reading 



\7o 



Here the time depending slope at y = yo is given by 



1 

*=2 



JdtAi 
J dtA2, 



(B.2) 



(B.3) 



meanwhile, the curvature of N(y) driven by the diffusive 
(stochastic) term in Eq. ( IB. II ). irreversibly decreases in time after 



In 10 

4 Jdt A 2 



(B.4) 
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from the large initial values corresponding to the initially mono- 
energetic distribution. Correspondingly, the time-dependent 
height at y = y follows 



No 



1 



1 



2 V^To 



JjdtA 2 



exp 



[fdtAi + JdtA 2 f 



4jdtA 2 



(B.5) 



To wit, Eq. (1B.2I > describes the evolution of the electron 
distribution, growing broader and broader under the effect of 
stochastic acceleration, while its peak moves from yrj to the cur- 
rent position 



7 m = 7o 



Jdt(M-h_) 



(B.6) 



under the contrasting actions of the systematic and stochastic 
accelerations. An important quantity to focus on for the emission 
properties is the r.m.s. energy 



Yp = 



j y 2 N(y) dy 
A| fN(y)dy 



7o 



J (A l+ 3A 2 )di 



7o 10* 



7m 



which is also the position for the peak of the distribution y 2 N(y). 

We have already derived in Appendix [A] the shapes of the 
spectra (synchrotron, and IC in both the Thomson and KN 
regimes) emitted by the distribution given in Eq. (01; here we 
stress the time dependence of their main spectral features. We 
can write for the synchrotron emissiorQ 



S oc 



„2 f (A l+ 3A 2 )dt 



y\ Vr, £, oc e 



2 f(Ai+5A 2 )dl 



A 2 dt 

for IC emission we have in the Thomson regime 

g 4 j(A ]+ 3A 2 )dr 



ccy2l07;(B.8) 



C oc 



4 / (A i +5A 1 )dt K yl _ 



JJ A 2 dt 

and in the extreme KN regime 

1 r i 

C oc — == oc Vr, £ oc eJ oc y p 10" 



(B.9) 



(B.10) 



dt 



Note that during flares, since A\ » A 2 , we expect the curvature to 
vary so little that we can therefore approximatively write S oc £ 
(see main text). 

Now we focus the above relations for the simple case of time 
independent A\ and A 2 , when J A\ 2 dt « A\ 2 t\ then we have 



2 1 A 2 



In 10 0.58 



(B.ll) 



A Ait 

while for the r.m.s. energy we have 

y p =yoe^ +M * )t ; (B.12) 
so for the peak frequency and the spectral curvature we have 



lo g yH=20li+5A2)f, 



* s= iM 5+ £)io g (iy 



3 Note that singular behaviors of 5 and C for t = are a consequence 
of the (differential) definition of the SED as vF v , and relate to the initial 
singularity of the particle energy distribution. Integrated quantities like 
F = J dv F v behave regularly. 




18 20 
log v (Hz) 



Figure B.l. Example of time evolution of the y 2 N(y) distri- 
bution (with peak energy y p ) due to stochastic and system- 
atic accelerations (upper panel) with A\ — 5 1CT 3 days -1 , A 2 
510" 4 days -1 



from the initial value y p = 10 3 , and of the re- 
lated evolution for the SSC SEDs (lower panel). In terms of the 
stochastic acceleration time t 2 , the time interval between each 
pair of lines is t 2 — t\ = 10~ 2 t 2 , corresponding to an observed 
time interval of about two days. 



where is a normalization frequency. It is seen that soon after 
the injection the curvature b s (proportional to r) drops rapidly, 
then progressively decreases more and more gently, while the 
peak frequency still increases. 

The value of A 2 can be evaluated from observing the syn- 
chrotron spectral curvatures b 2 and b\ at two times t 2 and t\, 
respectively (recall that b s * r/5); denoting with t 2 - 1\ this time 
interval, we have 



A 2 = 



0.58 
h-h 



1 +z 



(B.14) 



on the other hand, form observing the related synchrotron peaks 
£2 and £1, the value of A\ can be evaluated as 



(B.13) h = 



1 



2(f 2 -fi) h-h 



1+z 



(B.15) 



For example, in the case of Mrk 501 in the states of 7 and 
16 April 1997, we obtain (on assuming B * const) A\ - 
(2.3 + 1.1) yr 1 and A 2 = (1.8 ± 1.7)10 -1 yr 1 , correspond- 
ing to acceleration times T\ - l/Ai = (4.3 + 2.0)10 _1 yr and 
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Figure B.2. Example of time evolution of the y 2 N(y) distri- 
bution (with peak energy y p ) due to stochastic and system- 
atic accelerations (upper panel) with A\ = 5 1CT 3 days -1 , A 2 = 
5 1CT 4 days -1 from the initial value y p = 5 10 4 , and of the re- 
lated evolution for the SSC SEDs (lower panel). In terms of the 
stochastic acceleration time t 2 , the time interval between each 
pair of lines is f 2 — t\ = 1CT 2 T2, corresponding to an observed 
time interval of about two days. 



t 2 = 1/A 2 = (5.5 + 5.0) yr. Note that with the current data the 
evaluations of Ai and A 2 turn out to be affected by uncertainties 
considerably larger than the single curvatures b\ = 0.161+0.007 
and b 2 = 0.148 + 0.005. 

If the total energy available to the jet is limited (e.g., by 
the BZ limit, see text) we expect that J dt A\ and J dt A 2 can- 
not grow indefinitely, but are to attain a limiting value. At low 
energies where J dtA\ » f dt A 2 holds, we have S oc f and 
C oc e as before; at higher energies when J dt A\ reaches its 
limit, J dtAi «: J dt A 2 holds, leading to S oc £°- 6 and C oc e 6 . 
Eventually also J dt A 2 reaches its limit, and both the fluxes and 
the peak frequencies cannot grow any more. 

We add that the effect of radiative cooling is to reshape the 
distribution by moving high energy particles to lower energies; a 
corresponding deformation takes place in the emitted spectra at 
their high energy end, leaving nearly unaffected the peak fluxes 
for long cooling times. Then a term 77 ^ (y 2 N^j is added on the 
r. h. s. of Eq. ( IB. Il l, where for synchrotron radiative cooling 77 = 



the Kaplan approximation (Kaplan 119561) when A\ and A 2 are 
time independent; the solution is a power series in the parameter 
<7 - i] Jo / (2 A 2 ) <c 1, that reads 



N(y,t)*N ij-J 
where 

<r> = ro 10^ 

is particle mean energy. 



(B.16) 
(B.17) 



Appendix C: Pref actors for Eqs. © - (TT3) 

The HSZ SSC model is, as stated before, characterized by five 
parameters: the r.m.s. particle energy y p , the particle density n, 
the magnetic field B, the size of the emitting region R and the 
beaming factor 6. So the model may be constrained by five ob- 
servables that we denote with 

S C 

$i = : pro ' Cj = ; pro ' (C.I) 



10' 



erg 

cm 2 s 



10 J 



erg 

cm 2 s 



6 = 



£h = 



(C2) 



10* Hz ' 10 A Hz ' 

where the indexes i, j, k, h express the normalizations as demon- 
strated below; in addition, we denote with Atd the time in days 
for the source variations, and with D the distance of the source 
in Gpc. 

In the Thomson regime we find 
B = 10 -1 x 



13.5 x 

[2)5 b~i (1 +z)5 e 22 5 S-u^u' 1 C^rUr^] 



(C3) 



(C.4) 



■^cttcB 2 B 2 holds. The equation can still be solved on using 



R = 3.5 10 16 x (C.5) 

x [dH"5 (1 +z)- 1 2 en*S-u^ Af rf ^i 4 _1 C- U ~*] cm 

n = 5.7 x (C.6) 

x [b* D* (1 + j)5 ltf*" 1 ) ^ 4 2 C-n* S-n-i Af^-i] cm -3 

y p = 2.7 10 3 [loK 1 "*) e 22 i f 14 -i] . (C.7) 

In the extreme KN regime we obtain 

B = 1.7 10~ 2 x (C.8) 
x [2)Hi^l + z )-?£ 18 §S-io^ e26~^ C-ir^A?^] G 

5 = 4.2 x (C.9) 

x [dH-to lotC*- 1 ) (1 + z )i e 2 ^ S -wUi&^ C-u^ At d -i] 

R = 1.1 10 16 x (CIO) 
x [2)5 b~T° lO^i -1 ) (1 + z)~? e 2 6^ S-w* Af rf 5 f lg - l C-iT^J C m 

n = 9.5 10 2 x (C.ll) 

Tl 4 13/i 1 \ 4 LI 7 4 9 1 1 -1 

X [^3 z)-5 lO-"- 1 -*)^ + z )5£ 18 -c_ n 5 e 26 -5S_ 1() -5 Af d "5j cm" 3 

y p = 6.1 10 5 x (C.12) 
x [fera D~3 loM 1 -^) (1 +z)5 f 18 5 e 26 l c_n' A^5 S_ 10 _ §] . 
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Appendix D: KN correlations in the decay stage 

Here we briefly comment on the decay phase of a flare, when 
fluxes and frequencies decrease; the aim is at understanding 
the quadratic correlations between the y-ray and the X-ray 
fluxes observed in HBL objects like Mrk421 and Mrk501 (see 
Katarzynski et al. 2005; Fossati et al. I2008I ). Electron energy 
losses by adiabatic expansion (at a constant rate) can be ac- 
counted for by adding to the Fokker-Planck equation a nega- 
tive term that combines with the systematic acceleration rate 
(Kardashev 1 1 9621 ) into an overall coefficient At — > A\ — j. 

Correspondingly, the size of the source grows as R oc t\ 
with the total particle number remaining constant, in the extreme 
Klein-Nishina regime Eqs. lB.8l and lB.10l mav be rewritten as 



on neglecting the stochastic acceleration contribution. 

If both front shock thickness and gas clouds size scale as R( f) 
then A\ = 1 /(A t) holds. With A > 1 both the synchrotron and 
the IC flux will decrease. The particular choice A — 2 yields a 
quadratic correlation between the IC and the synchrotron peaks, 
that is, A log C — 2 A log S . One should not quail at the fine 
tuning corresponding to A = 2 literally, because lower values of 
1 < A < 2 are more effective in a mild KN regime. On the other 
hand, A cannot be large lest the shock is stalled. 




(D.l) 




(D.2) 



